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Abstract. We describe the ground state of a large, dilute, neutral atom Bose- 
Einstein condensate (BEC) doped with N strongly coupled mutually indistinguishable, 
bosonic neutral atoms (referred to as "impurity" ) in the polaron regime where the BEC 
density response to the impurity atoms remains significantly smaller than the average 
density of the surrounding BEC. Wc find that N impurity atoms with N 1 can 
self-localize at a lower value of the impurity-boson interaction strength than a single 
impurity atom. When the 'bare' short-range impurity-impurity repulsion does not play 
a significant role, the self-localization of multiple bosonic impurity atoms into the same 
single particle orbital (which we call co-self-localization) is the nucleation process of 
the phase separation transition. When the short-range impurity-impurity repulsion 
successfully competes with co-sclf-localization, the system may form a stable liquid of 
self-localized single impurity polarons. 
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1. Introduction 

Cold atom traps offer intriguing advantages for studying strong interaction physics in 
quantum many-body systems |1]. We focus on tlie cold atom prospect for conducting 
strongly coupled polaron studies. The polarons we consider are neutral atoms (referred 
to as "impurity atoms" in this paper) distinguishable from the BEC atoms, i.e., atoms 
in a different spin state or atoms of a different species embedded in a large dilute gas 
Bose-Einstein condensate (BEC). The BEC polaron is the impurity atom accompanied 
by its phonon cloud created by the interaction of the impurity with the BEC atoms. 
By varying an external magnetic field near a specific Feshbach resonant value [21 [3l H] 
experimentalists can tune the strength of the impurity-boson interactions and bring the 
polaron into the strongly coupled regime. In that regime, the BEC-polaron self-localizes 
[SI E] : the ground state properties become consistent with an impurity wavefunction of 
finite extent. The recently demonstrated species-specific potentials (e.g. the impurity 
atoms), such as laser generated potentials experienced by one kind of atoms [TJ [8], 
provides unparalleled prospects for the manipulation of such polaron objects. 

Attractive impurity-boson interactions give a collapsed ground state in 2D and 3D 
[9]. The cold atom realization of self- localized polarons with repulsive impurity-boson 
interactions can demonstrate polaron physics in regimes that are difficult to access in 
condensed matter [10]. Such cold atom experiments can then address fundamental 
questions pertaining to the nature of self-localization, the description of heavy polarons 
with a mass that is comparable to the energy scale of the boson modes and the many- 
body structure of systems with many polarons. At sufficiently low temperatures and 
sufficiently high polaron densities, cold atom traps may realize polaron superfluids 
without the complications that arise in describing the inevitable Coulomb interactions 
in electronic polarons 

Arguably, the polaron system provides the paradigm of a strongly interacting many- 
body environment modifying the properties of a single particle [10]. In the original 
polaron problem - the description of an electron in an ionic crystal - Landau [12] 
and later. Landau and Pekar [13] showed that the electron self-localizes in the strong 
coupling regime. In an attempt to understand the role of backflow in the roton feature 
of the helium superfluid dispersion, Feynman and Cohen studied a different polaron: 
an impurity atom in a condensed ^He superfluid [H]. The paper by Pines, Miller 
and Nozieres [I5] explicitly referred to this system as a polaron and pointed out the 
similarity between Feynman's ansatz wavefunction and the wavefunction implied by a 
second-order perturbation calculation. In this sense, the polarons here can be defined 
as the excitations created through the interaction of the impurity with the BEC atoms. 
Cold atom physics is now poised to realize polarons of this type at strong coupling, 
creating structures that are smaller than the coherence length of the BEC. These may 
provide a new handle with sub-coherence length resolution. 

The description of impurity-BEC polarons brings in aspects of related problems. 
Self-localization occurs at intermediate coupling, which is the challenging regime to 
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describe. A Feynman path integral description (TB] of the BEC-impurity polaron showed 
a surprising, non-monotonic dependence of the polaron size with coupling strength near 
the coupling strength at which the strong coupling description predicts self-localization 
[17]. In this paper, we will work within the strong coupling description in the assumption 
that this description correctly captures the important trends. We describe the co-self- 
localization of multiple impurity atoms in which the many-body ground state properties 
become consistent with bosonic impurity atoms occupying the same single particle 
wavefunction of finite extent. The starting point of and p], the assumption of a 
product state, was designed for the description of electron bubbles in helium superfluids 
|18] . Below, we show that the linearization of the BEC response to the impurity 
wavefunction reproduces the results of the Landau-Pekar elimination of the boson 
modes which couple as in the Frohlich Hamiltonian to the single impurity atom in 
the Bogoliubov approximation of the surrounding BEC. That description shows that 
the energy reduction produced by the BEC density response can be cast into the form 
of an attractive self-interaction potential. The linearization also shows that the Frohlich 
impurity-BEC coupling, the Landau-Pekar elimination and the resulting self-mediated 
interaction fail when the BEC-density response becomes comparable to the density of 
the surrounding BEC. When the BEC-density is not affected much by the localized 
impurity, the structure of the impurity-BEC system is that of a polaron. When the 
BEC density is expelled by the impurity, the impurity becomes a 'bubble' separated 
spatially from the BEC. The physics of multiple, bosonic self-localized impurities is 
connected with the phenomenon of zero-temperature Bose-Einstein condensate phase 
separation [19], [20], [21], as we point out below. The response of a dilute BEC to 
a localized object was studied in the context of a moving object that is significantly 
smaller than the healing length of the BEC p2] and in the context of a trapped ion 
interacting with the BEC atoms [23] via the polarization interaction of charged and 
polarizable particles. The polarization potential can also give rise to self-localization of 
the ion as shown in [23]. 

In this paper, we describe the ground state of a large BEC with multiple, strongly 
coupled, mutually indistinguishable, bosonic impurity atoms embedded in a large dilute 
gas Bose-Einstein condensate. In the absence of the impurity atoms we take the BEC to 
be homogeneous. The boson-boson, boson-impurity and impurity-impurity interactions 
are described by repulsive short-range effective interactions (contact interactions). 
As mentioned, we assume that the impurity atoms co-self-localize: the self-localized 
impurity atoms do not only overlap, they occupy the same single-particle orbital. We 
find that multiple impurity atoms can self-localize at smaller values of the impurity- 
boson coupling constant, than a single impurity atom can. We identify the minimal 
number of impurity atoms needed to stabilize the self-localized polaron structure for 
fixed coupling constant as a property of the droplet formation process in the nucleation 
of the BEC phase separation transition. When the short-range impurity-impurity 
interaction is sufficiently large, its competition with the self-localizing BEC-mediated 
interaction can lead to a regime in which one impurity atom self-localizes but not two 
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or more. 

The paper is organized as follows: in section [2l we discuss the magnitudes, 
dimensions and coupling constants of single BEC-impurity self-localization. In section |3l 
we use the Landau-Pekar method for describing the single neutral atom BEC-impurity 
problem so that we can later proceed beyond the homogeneous phonon description and 
investigate multiple impurities in the strong coupling limit. In section |U we examine 
the system of strongly coupled multi-impurity boson particles, followed by section |5l 
where we examine the multi-impurity polaron structure variationally and discuss the 
connection with BEC-phase separation dynamics. Finally, in section |6] we conclude. 

2. Single impurity self-localization in BEC 

First, we consider the zero-temperature self- localization of a neutral impurity atom 
embedded in a dilute gas BEC. The impurity atom couples to the elementary BEC 
excitations and forms a polaron. Strictly speaking, the resulting object resembles a 
piezo-polaron since the BEC excitations are acoustic phonon modes. The term 'piezo- 
polaron' was originally used by Mahan and Hopfield [25], [26], for conduction electrons 
interacting with the acoustic modes in a piezoelectric crystal [27], [28]. To bring the 
impurity atom to the strong coupling regime of self-localization, we assume that either 
the impurity-boson scattering length qib is Feshbach tuned to a large positive value 
or the boson-boson scattering length gbb is tuned to a small positive value. We 
describe the orders of magnitude and discuss the experimental challenge of avoiding 
large BEC depletion when impurity self-localization occurs. Since many BECs can 
be described as locally homogeneous (i.e., they are in the Thomas- Fermi limit), we 
consider the ground state of a homogeneous BEC of Nb indistinguishable boson atoms 
of mass niB, confined to a macroscopic volume Q, corresponding to an average density 
Po = Nb/^- Two bosons at position Xj and Xj effectively interact via a contact potential 
VBsi'^i — Xj) = [Airh'^aBB/'m'B] S (xj — Xj). The BEC is dilute, so that the gas parameter 
7 — V PoO-BB is niuch smaller than unity, i.e., 7 <^ 1. The chemical potential of the 
BEC bosons is /x^ = XbbPo and the BECs healing length - the distance over which 
the BEC-density tends to its asymptotic value if the condensate vanishes at a planar 
boundary - follows from equating the corresponding kinetic energy to fiB- We will work 
with a coherence length ^ for which /[2'mB^^] = 2fiB, which gives ^ = y/lGirpoaBB- 
The excitation of a collective phonon mode of momentum q then costs energy 

huj^ = h\q\c^l + {\q\0' (1) 

where c = \J p^b/^^b denotes the BEC sound velocity. 

Throughout the paper, we refer to the BEC atoms as "B" and to the impurity 
neutral atoms as "I" . When impurity atoms of mass m/ are added to the BEC, 
the impurity atom at interacts with a BEC atom at as well as with another 
impurity atom at via short-range interaction potentials Vib (fj — r^) = A/^^ {^i ~ ^j) 
and vii{Yi — Yi) = A//(5 (rj — r;), where Xjb = 2nh^aiB {m^^ + rrij^^ and A// = 
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[Anh^au /mi) and a/^, o^u are the scattering lengths describing the low energy impurity- 
boson and boson-boson scattering processes. We take all inter-particle interactions to 
be repulsive, i.e., aiB,ciBB,cin > 0, where qbb is the scattering length of the potential 
^BB = {4:7rh'^aBB/fnB)- If the impurities are bosonic and coexist with the BEC-bosons 
at sufficient density and at sufficiently low temperature, the system is a mixture of a 
"B" and an "I" condensate. The "B"-BEC mediates interactions between the impurity 
atoms that are described by an attractive Yukawa potential of range 

Vs{r, - r,) = -—^e-d"--^^!/?) (2) 

where we characterize the strength by an effective charge Q to emphasize the useful 
analogy with Coulomb interactions. The appearance of a Yukawa interaction is not 
surprising. It is known in quantum field theory [21] that the interactions mediated by 
a scalar boson field (provided by the BEC) takes the form of an attractive Yukawa 
potential in a non-relativistic case. In accordance, the impurity atoms in the "I"- 
BEC distributed at density p/ (r) experience a local chemical potential /i/ (r) with a 
Pi dependent contribution 

f^i (r) = XiBPo + hiPi (r) - j dh' [QV |r - r'|] e-l-^'l/^/ (r') 

= XibPo + (A// - 4<2g2^) 

if the impurity BEC is distributed homogeneously. 

The above chemical potential term implies a diverging "I" compressibility (~ 
[dfii/dpi]"^) when A// — Ati^'^Q^. It is the diverging compressibility that causes the 
zero-temperature BEC phase separation that we describe in section |5l This transition 
occurs when the Fetter-Colson relation [19] is satisfied Xn — Xib {Xib/ Xbb), where 

XbB = {4:'7lh'^aBB/TTT'B) SO that 

XlB 1 . / XiB 



Q' = ^IBT^J-r2 = ^^^aBB (4) 

Xbb^t^K \XbbJ 
is the value of the Yukawa strength parameter. 

The BEC mediated interaction is a consequence of BEC density correlations: an 
impurity atom i at causes a density variation in the BEC that is experienced by 
impurity atom j at rj as a change in its BEC mean field energy. Describing the impurity 
atoms by a quantum wavefunction, we find that even a single impurity is affected by 
BEC density correlations. The impurity interaction energy density at r is influenced by 
the impurity density at r' as the BEC responds to the density of the entire impurity 
wavefunction. As we will see below, a single impurity then experiences an effective 
self-interaction Vs (r — r') = — Q^e"''"^'"'!/^/ |r — r'|. The integration of Vs (r — r') over 
r and r' weighted by ^pi{r)pi{r') gives the gain in impurity-boson interaction energy 
by the adjustment of the BEC density to the impurity wavefunction. 

In the strong coupling limit, the effective self-interaction exceeds the kinetic energy 
cost of localizing the impurity. For the Yukawa self-interaction to overcome the kinetic 
energy, the impurity has to localize to a size comparable to or less than ^. In that case. 
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the self-interaction is Coulomb-hke in most of the impurity-region and, hence, is efficient 
at binding. As in the Hydrogen atom description, the length scale Rq is determined by 
comparing Coulomb and kinetic energies 

Ro niiRl ^ ' 

Using equation (jl]), /i^ = ^^/[4mB^^], and = IGnpQaBB, we find that the effective 
Rydberg length is equal to 

^11 ^ 

° 47rpoa?ij f 1 I m.\ (i + ms) 

\ ruB J \ mi J 

Aside from the mass factor (^1 + ~^ ™?^) ' mean free path of 

an impurity atom that encounters hard-sphere scatterers of radius a/s distributed 
at average density po- The energy scale of impurity self- localization is given by the 
corresponding Rydberg energy 

^TT^ = 2 Urttpd 1 H 1 H 



The mediated self-interaction can only bind with sufficient efficiency if the natural 
impurity size Rq is comparable to or shorter than ^. The ratio of the BEC-coherence 
and Rydberg lengths, C,/Ro, is the only parameter left in the description of the impurity 
wavefunction after scaling the energy and length by Eq and Rq, respectively. We call 
this dimensionless ratio, the impurity-boson coupling parameter, 

^.X.V^.Sfl + ^VlH-^V (8) 



Ro y asB \ rriB J \ mj 

Calculations [5] indicate that a single impurity atom can self- localize if /3 > /3c (1) = 4.7, 
where /3c (1) stands for the critical /3 for a single impurity atom. As /3 further increases, 
the self-localized impurity size shrinks to become a point-like object (smaller than the 
BEC coherence length). The calculation of a point-like potential moving with velocity 
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V [22] gave an energy m£)f^/2 with m£, = \raB\fTi\j , so that we might expect the 
effective polaron mass m* to be equal to m* = m/ -|- m^) or 

m* 2/3 

— = 1 + — ^ (9) 



However, the calculation in [22] was based on perturbation theory so that its validity 
does not extend to the large coupling limit. For the classic polaron treatment that 
describes an electron coupled to optical phonons, Feynman estimated an effective mass 
that varies as the fourth power of the coupling constant [TO] . 

It is useful to express the observables in terms of the impurity-boson coupling 
constant /3. For instance, 

g2 ^ 2^/3Vi?i?o (10) 
mi 

Eo = 2^/3Vb (11) 



Multi-impurity polarons in a dilute Bose-Einstein condensate 



7 



with To = h/Eo representing the time scale relevant to impurity self-localization. 

The polaron description of the BEC-impurity assumes that the phonon modes 
are excitations of a homogeneous BEC. However, as the impurity-boson interaction 
is increased to reach the critical coupling f3c{l), the impurity can deplete the "B" 
condensate in its vicinity. The above description breaks down when the local 
BEC-density variation 6pB becomes comparable to pQ. For repulsive impurity-boson 
interactions we expect the single impurity wavefunction to take on a phase-separated 
bubble profile when 6pB — Po- While this is interesting by itself, the impurity loses its 
Frohlich-coupled polaron character in the bubble limit. As the self- interaction energy 
is the change in boson-impurity energy caused by the impurity induced variation in 
BEC-density, Xib^Pb ~ Eq or 



5pB Eo EoXbb 4^1/4 /^^l^i (12) 



Po A/spo Pb Xib V 



where 7 represents the BEC gas parameter, 7 = a/ Po<^bb^ defined earlier in this section 
(above equation ([1])). 

Can the self-localization of a single BEC-impurity be realized in cold atom 
experiments? At the end of the next section, we discuss whether cold atom technology 
can access the polaron regime - the parameter region where the density variation, 6pB, 
in equation (fT2l) remains significantly smaller than pQ. Here we discuss the experimental 
accessibility of the strong impurity-boson coupling regime, /3 > 5. Three experimental 
challenges must be met. Two of them relate to the requisite increase in impurity- 
boson scattering length. First of all, if that increase is effected by a magnetically 
controlled Feshbach resonance, the magnetic field has to be sufficiently homogeneous 
and constant to avoid significant variations of the impurity-boson interactions. Secondly, 
three body recombination into the large, two atom impurity-boson bound (dimer) state 
that becomes degenerate with the incident channel impurity-boson continuum as ajB 
diverges, should not occur before the BEC impurity has self-localized. As a third 
condition, we mention that the BEC-temperature should be sufficiently low to prevent 
thermal fluctuations from de- self-localizing the strongly coupled BEC impurity. 

How large does ajB have to be? Estimating the necessary ajB from 



ajB 

O-BB 



/3/7 



\ ^ V rriB J \ mi 



(13) 



with 7 ~ 10 ^ and a realistic mass ratio of ~ 10, we find that ~ 10^ if /3 ~ 10. The 

' ' ass 

Feshbach scattering length with magnetic field dependence is given by 

aiB = ajB,o 1 - p , (14) 

B — Bq_ 

where 07^,0 denotes the background scattering length, A is the resonance width, and 
Bq is the resonant field strength. Hence, a slow magnetic field variation 5B induces an 
aiB variation 

A 

5aiB = o-iBfi— r—25B. (15) 

[B - Bo) 



Multi-impurity polarons in a dilute Bose-Einstein condensate 



8 



For B near Bq, a^s — — fl/_B,oA/ {B — Bq), so that 



6a J B ajB 



(^IB CllBfi 



6B 



A 



(16) 



Requiring \SajB/ajB\ < e leads to \SB/A\ < e {aiB,o/(^iB)- Assuming, for instance, 
aiBfi/aiB = 1/30, e ~ 0.03, we obtain \SB/A\ < IQ-^. For A = IG, the magnetic field 
variation is 6B < ImG, which can be achieved in today's laboratories. In contrast, a 
broader resonance may require SB < lOmG (A = lOG) or SB < O.IG (A = lOOG). 

We now estimate the recombination time of the BEC impurity atom in the 
worst case scenario. Applying the same reasoning as in [30] where the three-body 
recombination rate for a BEC of large scattering length was first calculated, we expect 
a recombination rate of the form 

- - '4Jl + 2(l + ^)±aWo. (17) 
Tr y/3\l V mi J TUB 

For a BEC of large and positive a^B, substituting mi ttib and ajB ^ clbb equation 
(fT7|) gives an overestimate of the recombination rate (see [31]). The authors of reference 
[3T| found that the power has to be multiplied by an oscillating function of magnitude 
less than unity. Near a node of the oscillating function, the recombination rate is much 
reduced. Nevertheless, using equation f lT7|) . we find 

^« = 'W / (18) 

V 

where tb = h/ ^b represents the time scale on which the BEC can respond to the 
impurity and form the phonon cloud that accompanies the impurity in the polaron 
state. Note that for mj ~ lOm^, /3 ~ 10, tr can exceed tb by more than two orders 
of magnitude. Finally, we comment on the temperature requirement for observing BEC 
impurity self-localization. Determining the polaron properties in a finite temperature 
Feynman path integral calculation, reference [16] identified polaron self-localization with 
a sudden and unexpected increase in polaron size as /3 increases near the critical value 
obtained from the strong coupling description. The polaron size increase was observed 
for a ^Li-impurity in a ^^A^a-BEC in the case where the temperature (/^bT) was smaller 
than half the boson chemical potential. Therefore, we expect that self-localization takes 
place when /c^T < yU^/, where / can be as large as 0.5. For a BEC of = ^ x 300iJz, 
the temperature requirement of T < 7.5nK can be easily achieved in today's cold atom 
experiments. 



3. Landau-Pekar elimination of phonon modes 

References [5] and [6] describe the self-localized BEC-impurity state as a product 
wavefunction. Specifically, the many-body state of an impurity of position r and Nb 
bosons of position xi, X2, . . . xat^, Bose-condensed in volume Q so that the average boson 
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density is po = Nb/^, was represented by a wavefunction ip (r) x (^i) X2 (X2) • • • X (^iv)) 
where ^^id x denote the normahzed single-particle impurity and BEC boson 
wavefunctions, respectively. In the absence of an impurity, the BEC-field (r) = 
y/NsX (r) is V9 = s/po, position independent for a homogeneous BEC. Assuming that the 
BEC response to the impurity remains smaller than ^/po, (p = ^/po + Sip, the equations 
can be linearized in 6ip. The resulting description is equivalent to the Landau-Pekar 
treatment we describe below. The product wavefunction reveals the connection with the 
analogous "fluid with embedded droplet" system. The linearization procedure indicates 
how to proceed beyond the homogeneous phonon description. On the other hand, the 
Landau-Pekar method is most easily generalized to describe multiple impurities in the 
strong coupling limit. 

In the absence of an impurity, the excitations of the dilute BEC are well-described 
in the Bogoliubov approximation. This description transforms the operators that create 
and annihilate boson particles of momentum q, 6^, 6q to quasi-particle (phonon) creation 
and annihilation operators a|j, Oq, that diagonalize the linearized Hamiltonian. Adding 
a constant term, the transformed BEC-Hamiltonian reads 

^/kUqUtaq + i) , (19) 



q 



B- 



where hwq is stated in equation ([T]). The first step of the Bogoliubov procedure 
replaces the fcq=o; ^q=o-operators by V^iVg and expands the energy in powers of \/Nj 
Accordingly, the boson density operator becomes 

PB,q ^ A^B^q^o + 




= Nb5,,o + y^f-^-^^^ («i + (20) 

where the last step followed from implementing the Bogoliubov transformation. 

Introducing the impurity density operator p/^q, the impurity-boson interaction takes 
the form 



7^ 2^P/,-qPB,q - ^IBpOni+ ^ 2^ y ^ P/,_q (fl^ + a_q) , (21) 

where fij represents the number operator of the I-atoms, p/,k=o = ^i- The second term 
on the right hand side resembles the electron-phonon interaction term in the Frohlich 
Hamiltonian 

^ '^-q^^-q (^q + ^-q) (22) 
^ q 

where Mq = Xisy/Po and i^_q = v\q\ = ^/KhPg"^ /2mB) /hujq\- The phonon part of the 
impurity-BEC Hamiltonian, Hp, reads 



^t-^ ^ 1\ y2Mo . M ^ 
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The strong coupling approximation assumes that the impurity is described by a well- 
defined wavefunction that is not entangled with the boson B state. As a consequence, 
the coupling to other impurity states can be neglected and the impurity density operator 
P/^q can be replaced by its expectation value, p/^q — )■ p/^q = / d^r exp(— iq-r) {ip {r)f. To 
bring out the oscillator nature of the boson modes, we introduce the oscillator coordinate 
and momentum operators that are equivalent to the creation and annihilation operators, 

^'^=^^ 

at — d^a 

^ (23) 

In this notation, the phonon Hamiltonian Hp becomes a sum over second order 
polynomials in 0, 



q 



5q<Pq + TT^TTq j + — ^I^„qP/,„q</)q 



(24) 



Completing the squares in equation (^^, we write the phonon Hamiltonian as a sum 
over displaced oscillator Hamiltonians 



q q ^ 



where the ^-operator denote the displaced boson oscillator coordinates 

As the displaced ground state energy is the same as that of the original oscillator, the 
energy reduction caused by the phonon response to the impurity density is given by the 
remainder 



M2 



q ^ 

Using p/^q = f d'^r exp [zq ■ r] |^/;(r)|2 this expression can be cast in the form of an 
interaction energy 

AE, = i y dh j rf^r'l^ (r) I V, (r - r') |^ (r') \\ (28) 

where Vs represents the inverse Fourier transform of —2MQi'^/hL0q. The expression 
holds regardless of the momentum dependence of Wq and z/q which are different for ionic 
crystal (optical) and piezoelectric (acoustic) polarons. In the case of the BEC-impurity 
or "Bogoliubov polaron" , Vg is the attractive Yukawa interaction, 

1 f /-2MoV2 



{2Tiy J ' \ hw^ 

\ -^/B 1 1 f ,3 (-iq r) \ 



(29) 
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where we used equation (jl]) to obtain the above expression. The total change in ground 
state energy E caused by doping the BEC with a strongly coupled impurity atom, 
includes the kinetic energy cost of localizing the impurity is 



E = XiBPo / d'r |^(r)|^+ / d'r r {r 



d-^v / d^v'\i){r 



2m/ 



(30) 



The optimal wavefunction ip follows from minimizing E while ensuring that ip (r) 
is normalized. We include the normalization constraint by introducing a Lagrange 
multiplier /i. The extremum of the associated functional F = E' — /i [J d^r (r)|^ — l] 
is found by requiring 6E/6ip* (r) = 0, which leads to 

f^ii^ (r) = V (r) 



+ 



2m/ 
dh' 



|r-r'|/C 



l^(r )l^^(r) 



(31) 



where fij = fi — XibPo- The impurity wavefunction that solves equation fl?T]) is localized 
if and only if fij is negative. Multiplying equation (!3T!) by ip* (r), integrating over r, 
and making use of the normalization condition, we obtain 



Pi 



+ / d^r / d:h' \^ (r 



2m/ 

|2 



^(r) 



\r~r'\/( 



|^(r')r 



(32) 



which differs from equation ( l30i) by the factor of 2 in the self-interaction term. 

The numerical solution of equation (!3T]) . as calculated iteratively in [5], shows 
that the profile of ijj generally resembles a Gaussian function. Assuming a normalized 
Gaussian wavefunction, ip [r) = e"'""' Z^*^^/ (ttct^)^^^, we solve the problem variationally 
by minimizing E with respect to a. To know if the variational description predicts self- 
localization, we test if at the minimal width a, pj (a) < 0. Working out the integrals 
we find 

~ 3 



where /, 



Eia) 



f{y) 



2a 



f 



duu exp 



2/exp 



u 



4 m/(T^ ' 



exp {-uy) 



(33) 



21 



erfc 



y 



(34) 



accounts for the reduced efficiency of the Yukawa self-interaction when the Gaussian 
width exceeds the BEC-coherence length (as f{y) — )■ when y ^ 1). If we scale the 
energies by the Rydberg energy of equation ((7]), — (E/Eq), pi (pj/Eq), and the 
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length by the Rydberg length of equation x = a/ Rq, /3 = ^/ Rq, then E and /i/ take 
the form 

'''' (35) 



fiiix) f^ffA 2^ 3 





Figure 2. The variationally determined chemical impurity potential fii. In the strong 
coupling description, the BEC-impurity polaron self-localizes when /i/ < 0. 

In figure [H we show the variationally determined Gaussian width a in units of the 
Rydberg length scale Rq. The extent of the impurity wavefunction is proportional to 
this width with a factor of proportionality that is -\/3/2, y/ < r'^ > = 3 / 2a. In the 
Gaussian variational description, the bound state minimum of equation fl35|l disappears 
when /3 < 5. In figure we show the variationally obtained fij. As fij becomes negative 
when (3 > 5.4, the Gaussian variational description predicts that the BEG impurity 
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polaron self-localizes when (3 > (3^ with (3^ = 5.4. Since /3c gives the minimal coupling 
for self-localizing a single impurity atom, = 1, we will denote as /3c(l) = 5.4. If 
the extent of the polaron is significantly smaller than the BEC coherence length cr ^ ^, 
equivalent to x <C /3, we can Taylor expand / as /(x//3) f» 1 — ^Jti /2{x/ (3). In this limit, 
the impurity system energy E/Eq = —^2Itxx~^ +3/ \lx^\ + 1/13 attains a minimum 
at X = 3-\/7r/2 corresponding to \/ < r"^ > = (3/2) \/37ri?o or x — )■ 3^/7/2 ^ 3.76, the 
large coupling limit that the graph of figure [1] tends to as /3 — ?■ oo. Within the same 
Taylor-expansion, fii/Eo = — y^277r(2/x) + (3/ [2x^]) + 2//3, takes on a negative value 
if 13 exceeds 27r. Therefore, in the small polaron size approximation, /3c(l) ~ Stt, as 
compared to the variational result /3c(l) = 5.4 and /3c(l) = 4.7, numerically determined 
by solving the non-linear equation, equation [3TJ 

The BEC density response to the self-localized impurity is a consequence of the 
oscillator (0) displacement in the Landau-Pekar description. We calculate the density 
variation in the Bogoliubov approximation. 



6pB (x) 



PB (Xj - po 



k 



-pi,k 



(37) 



Using ( $k / = we find 
6pB (x) 



Po 



Pol_ 
Mo ft 



ikx 



2M|z^ 



(38) 



Recognizing the expression in the square brackets as the Fourier transform of Vs and 
replacing J^jM^ = 1/\ib, we obtain 



6pB (x) 



dh- 



-|x-r|/? 



(39) 



XiB J |x - r| 

which is identical to the expression found in |5] by linearizing the BEC-response. Note 
that the Landau-Pekar description and the linearization of the product state give the 
same results for all observables. 

We now determine the peak value of the BEC density variation. If we choose the 
origin at the center-of-mass of the impurity in [^/^(r)!^, the peak density is reached at 
X = 0. In the true ground state, the impurity is not localized in any specific localization 
but the strong-coupling product wavefunction breaks translational symmetry. While the 
wavefunction does not have the correct symmetry, we know that some of the relevant 
polaron properties such as energy, mass, and polaron size give the correct values in the 
strong-coupling limit as they do in the traditional optical polaron problem. Assuming 
the width to take on the strong coupling value, a = 3^y7rj2Ro, 



SpB (x = 0) 
Po 



Eq Xbb 4^2 
Pb ^ib 37r 



/ 



2/3 J 



(40) 
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where the factor ^^ff- is given by equation (I12p . Assuming 3-^71/2/3 ^ 1, we find 



Spb (x = 0) 



16V2 



V nil V 



po 37r3/4 

In the polaron hmit, the hmit in which the Bogohubov- approximated impurity-boson 
couphng of the Frohhch type is vahd, the BEC-density response should be smaller than 
the average BEC density, Sps/po < e, where e is a fraction of unity. We find that the 
average boson BEC-density should not exceed a value pmax with 

1 (mBV 1 _ 0.96 X 

^"^■^^ alAmiJ (3.196)^/3^ a|« \mj J /je" ^^^^ 

Depending on the atoms trapped in the experiment, the above condition can pose a 
challenge. Taking the impurity atom to be a spin-fiipped atom in a ^'^Na-BEC, choosing 
e = 1/5 and assuming that ajB is Feshbach tuned while ass = CLxripiet = 3nm, we 
find Pmax = 3.6 X lO^'' cm~^, which is a rather low density to make observations with. 
On the other hand, embedding a ^^Rb-impurity atom into ''Li BEC of asB = O.Snm, 
with e = 1/5 leads to pmax = 1-2 x lO^^cm^ [/3//3c(l)]~^, yielding a maximal density of 
2 X lO^^cm"^ even for /3 = 10. 

4. Co-self-localization of bosonic BEC-impurities 

Most schemes for realizing strong coupling BEC- impurity polarons would create multiple 
impurities that are indistinguishable in the same cold atom trap. In this section, we 
describe the ground state of the impurity system in a homogeneous BEC. We assume 
that the impurity atoms are indistinguishable bosons and that the impurity-impurity 
interactions XjjS (r^ — rj) are repulsive. A// > 0. 

We will characterize the strength of the short-range impurity-impurity interactions 
by a third dimensionless coupling constant 

The system's behavior is determined by three dimensionless parameters: the BEC gas 
parameter 7 = ^/poo^, which affects the BEC-density response to the self-localized 
impurity but not the impurity wavefunction in the polaron limit, the impurity-boson 



coupling constant, /3 = ^/7l^y poaj^/aBB + rnB/mj) {1 + mj/rriB), and the above 
defined impurity-impurity coupling constant, which can also be written as 

ri = 8^2/37 (^] . (44) 
\aBB / 

We extend the above descriptions to treat A^ impurities of position rj, z = 1, 2, . . . A^, 
and Nb bosons of position Xj, i = 1,2,... Nb- The strong coupling ground state wave- 
function can be written as a product state "if (ri, r2, . . . fat) xb (xi) Xb {'^2) ■ ■ - Xb {'^Nb)- 
Alternatively, the Landau-Pekar treatment is easily generalized as the impurity-boson 
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coupling is described by the same Frohlich-like term of equation f l25p . In the strong- 
couphng hmit, the density operator is replaced by its expectation value Pq — )■ Pq where 

Pq = y dh, j dh2 ... j dh^ J2 ^''""^ 1^ (^1' • • • ^n)\' . (45) 

The strong coupling Landau-Pekar phonon elimination leads to a phonon-induced 
ground state energy reduction 

^E.p = ^j dhi j dh2 ...J dhr, J dh[ J d% ...J d'r'j, 

X 1^ (n, r2, . . . r;v) I' Yl (^^ - 1^ (^'i' • • • ^'n)\' ■ (46) 

We can argue that the i = j terms describe self-interaction while the i ^ j terms describe 
the BEC-mediated interactions. However, this distinction becomes less meaningful when 
the impurities are indistinguishable and \1/ satisfies the bosonic or fermionic permutation 
symmetry. In this case, all i,j contributions are equal and AEp becomes 



AEp = iiV^ J dh J dh'pr (r) K (r - r') p, (r') 



(47) 



The localization of the impurities comes at a kinetic energy cost and the impurity- 
impurity overlap comes at an interaction energy cost. We describe the direct impurity- 
impurity interactions by an effective potential vu (rj — r^) = A//(5 (rj — r^). Adding 
strongly coupled impurities to a homogeneous dilute BEC changes the ground state 
energy by an amount E where 



E = NXibPo 



+ J dhi J dh2 ... J d^N tp* (ri, r2, . . . rjv) ^ (^~) ^ ^^i' "^2' • • 
+ AEp 

+ [ dh, [ dh^... [ rf^r^v f dh[ [ d%... [ rfV^V^*(ri,r2,...r^) 



The mediated inter-particle interaction potential, K, is attractive [Vg < 0) and favors 
impurity overlap. Maximal overlap is achieved when the N bosons occupy the same 
single particle orbital, ip (^i, r2, . . . r^) = ip (ri) ip (r2) . . .'ip (fat). To emphasize that the 
boson impurity particles not only self-localize, but occupy the same orbital (localizing 
in Hilbert space), we will refer to the co-self-localization of the bosonic impurity atoms. 
In that case, p/ (r) = (r)|^ and using equation ( H7|) in the E expression we find 



E = NXjBPo + N I d'vr (r) ( ) (r) 



iV(A^-l) 
2 



\n fdh\i,{r)\ 



4 



j d^r j dh' \ij (r) I' K (r - r') |^ (r') |' , (49) 
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where ip (r) is assumed to be normalized. In determining the optimal '(/'-orbital, 
we introduce a Lagrangian multiplier Nfi to ensure J d^r\ip{r)f = 1. Minimizing 
E — Nji (J d^r \ip (r)l^ — l) with respect to ijj* yields 

/i,^ (r) = ^^^V^ (r) + (iV - 1) A,, 1^ (r) I' ^ (r) 

- NQ^ / d^r'f _ 1^ (r')l' ^ (r) , (50) 



where /i/ = fi — XjbPo- The initial A^-factor in the Lagrange multiplier ensures that 
fii becomes the impurity chemical potential in the limit of large impurity numbers A^. 
Solving the above equation for fij by multiplying both sides by ^Z^* (r) and integrating 
over r, we obtain 

d'rr (r) (^^] i' (r) 



+(iv-i)A,, y d\\i,{j)t 

-N J dh J rfV \4j (r')| ^ 1^ 1^ (r')| . (51) 

The physical fij value minimizes E in equation f H9|) . For negative values of fij, ijj 
becomes a localized wavefunction. At a distance exceeding ^ from where the impurity 
density becomes small, ip takes on the Yukawa profile of a bound state. 

5. Variational description of BEC-impurity co-self- localization 

Alternative to solving the non-linear Schrodinger equation of equation (l50ll . we can 
approximate the impurity orbital ip{r) variationally. As in section |3l we insert a 
normalized Gaussian wavefunction, ip{r) = e"^'^"^"^ (^^r^) into the energy E of 
equation (149|) and we minimize E with respect to the Gaussian width a. The treatment 
predicts self-localization of the A^ co- locating impurities if /i/ of equation ( I5T1) . evaluated 
at a that minimizes E, gives a negative value. 

5.1. The variational Junctionals 

In the scaling units, Eq, Rq, introduced above to treat the single impurity polaron 
problem, the E and /i/ expressions take the form 

and 



Eq X \ n \^ J 2x'^ x^ 

where x = c/Rq and where the impurity-impurity coupling was defined above (see 
equation (glD), r] = {XiiRq^) / [[2Trf^E, 
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Two interesting and relevant limits can be understood analytically: the 'pinprick' 
limit in which the size of the impurity wavefunction shrinks well below the coherence 
length of the BEC and the large impurity number limit in which the impurity 
wavefunction extent is determined by the competition of the BEC-mediated interaction 
and the impurity- impurity repulsion. 



5.2. The pinprick polaron limit 

The coherence length ^ is the natural length scale of the BEC. If an object has 
a subcoherence length size L, i.e., L <^ ^, the BEC response to its presence is 
indistinguishable from the BEC-response to a point object. In the presence of co- 
self-localizing boson impurities, the BEC-mediated impurity-impurity attractions that 
contribute a term ~ A^^ can further decrease the impurity size to be sub-coherence 
length — the impurity becomes a 'pinprick'. In this limit, the / argument, x//3 = a/^, 
is much less than 1 and / can be Taylor-expanded in the E and fij expressions, 

E (x; TV, /3, rj) = E, (x; TV, rj) + (54) 

and 

2NE 

fii{x;N,P,7]) = fi,{x;N,r]) + ^-^, (55) 

where Ec and Hc denote the Gaussian expectation values for Coulomb (instead of 
Yukawa) self-interactions, 

E{N) 12 ^8 N{N-1) 7] 



+ N 



Eq X \ n 2 x^ 

f^i (N) 2N [2 3 r/ 

In this limit, the equilibrium Gaussian width follows from the requirement dEc/dx = 
at fixed and 77, yielding x{N,ri). The critical coupling /3c {N,ri) follows from /i/ = 0; 

(3c{N,v)^ , , (57) 

He {x;N,r]) /Eq 

although this approach becomes somewhat questionable for determining the impurity 
properties near the critical coupling since the Taylor expansion is not generally valid 
near the self- localization point. Furthermore, if the short-range impurity- impurity 
interactions can be neglected {rj — ?■ 0) in the pinprick multi-impurity polaron, 
the condition, dE^/dx = leads to the equilibrium width x = "i^/n 12 /N ^ 

x{N = 1,?7 = 0) /N. The expectation value of Ec scales as iV'^ and the binding energy 
takes the form 

E {N- /3, 77 = 0) = -^Eo + ^E^- (58) 

OTT p 
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On the other hand, the expectation value of fic scales as A^^ , fic/Eo = —N'^/n and from 
equation f l57|) we estimate the critical coupling constant as 

/3c - |. (59) 

Hence, bosonic impurities can co-self-localize at much smaller values of the boson- 
impurity coupling strength than the critical coupling for single- impurity self-localization. 
The self-localized multi-impurity polaron can be significantly smaller than the BEC- 
coherence length. Specifically, casting the impurity wavefunction extent, \/Jr^ = 
1^3/2(7 = a/3 / 2xRo and the binding energy in terms of the BEC-length and energy 
scales, we find 



3V37re 

2nN (3' 

'niB 



For small, but finite impurity-impurity coupling 77, the impurity-impurity interaction 
becomes significant as increases. Solving for the equilibrium width x from 

dEjdx = 0, we obtain x = (3^/^^ (1/2) (^1 + ^Jl + (2/3) r]N{N - , so 

that the impurity-impurity repulsion significantly increases the extent of the impurity 
wavefunction when N{N — 1) > (^3y^7i/2/2r]j . In figure |3] we show the Gaussian 
width of the three- impurity polaron wavefunction at coupling /3 = 4 as a function 
of r], as determined numerically from the full functionals, equations (!52|) and ( !53|) . 
and as determined from the Taylor expanded approximation equations (1541) and ( l55l) . 
Notice that the Taylor expansion only works well for small values of rj as the impurity- 
impurity repulsion sufficiently increases the impurity radius to invalidate the Taylor 
expansion assumption ^ 1. As the expansion replaces the Yukawa self- interaction 
by a Coulomb attraction (more efficient at binding the polaron), the E^ approximation 
underestimates the size of the impurity wavefunction. 

The effect of the impurity-impurity repulsion is particularly striking in the 
dependence of the critical (3c on the impurity numbers A^, as shown in figure HI This 
figure shows f3c{N, rf) for A^ = 1 (independent of 77), 2, 3, 10, at 77 = 1, 10, 30 and 40. 
At ?7 = 0, [3c{N,ri = 0) decreases with increasing impurity numbers and l3c{N,'r] = 0) 
vanishes at large A^. At finite, fixed 77 with rj < 24 the critical coupling Pc{N, r]) (shown 
by dots in the same color) still decrease with increasing A^, but Pc{N, rf) tends to a finite 
value f3c{N — )• 00,7/) that depends on rj. At 7/ > 24, /3c(2,7/) > /3c(l), but then (3c{N,ri), 
N >2 still falls off monotonically with increasing impurity numbers A^. 

This dependence has interesting implications. If the BEC-impurity system has 
an impurity-boson coupling constant f3 in the range [/3c(A^^ — 1), /3c(^^)]5 the above 
description suggests that A^^ or more bosonic impurities would co-self-localize in the 
ground state, whereas Nf^ — 1 or less impurities would not. The minimal number of 
particles required to make a stable, self-localized structure is reminiscent of the familiar 
nucleation processes that take place in first-order transitions. In that case, the formation 
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Figure 3. The expected Gaussian width of the BEG impurity wavefunction that 
describes the ground state of the BEG with three bosonic impurity atoms occupying 
the same single particle (Gaussian) wavefunction, shown for /3 = 4 as a function of the 
impurity-impurity coupling constant rj. The upper curve, shown by the full line in blue, 
shows the variationally determined width, whereas the lower graph (shown in red dots) 
pictures the Gaussian width determined in the Taylor-expanded approximation. The 
Taylor expansion treats the BEG mediated interaction as a Goulomb interaction and 
underestimates the size of the impurity wavefunction, which grows as a consequence 
of the short-range impurity- impurity repulsion. 

of a stable droplet of the new phase also requires a minimal amount of matter in the new 
phase to overcome the relative cost of the surface energy (here, kinetic energy plays the 
role of surface energy). We tentatively identify the co-self-localization of impurities as 
the onset of the nucleation process of B-I phase separation, a connection we will further 
explore below. 
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Figure 4. Gritical values of the impurity-boson coupling constant /3, /3c, above which 
N bosonic impurity atoms self-localize in the BEG-impurity ground state described in 
the strong coupling approximation. The graph shows /3c for different impurity atom 
numbers N , N — 1, 2, . . . , 10 , set out on the horizontal axis, and gives the variationally 
determined /3c- values for different impurity-impurity coupling values, rj: red diamonds- 
77 = 0, black plus- 7] — 1, black square- = 10, black x- 77 = 30 and blue circle- rj — 40. 
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5.3. BEC-phase separation 

In a phase separation transition, initially miscible fluids become immiscible and separate 
from each other. At zero temperature, the ground state of a gas of I and B weakly 
interacting boson particles changes from a mixed to a separated configuration as 
the impurity-boson scattering length is increased. The natural coupling constant to 
characterize this zero temperature phase transition is given by 



As g varies from g < 1 (miscible mixture of I- and B- BECs) to g > 1 (ground state 
corresponding to separated BECs), the system undergoes a phase transition that appears 
to be second order - the compressibility and a coherence length diverge. If the I BEC 
is the smaller condensate, the ground state gathers the I-bosons into a spherical region 
as this geometry minimizes the surface energy. The surface tension itself vanishes as 
— 7- 1^, but is finite when g > 1. Hence, we expect separation dynamics at (? > 1 to be 
similar to that of a first-order transition. 

For large I- and B- BECs, the onset of the separation dynamics is described by 
the growing instability of the collective excitations of the homogeneous BEC-mixture. 
But how does the separation dynamics set in if only a few I-atoms are present and the 
I-bosons have not acquired macroscopic phase coherence, whereas the temperature is 
sufficiently low to treat the B-BEC as a zero temperature system? We suggest that 
when fluctuations gather a minimal number A'^ of I-bosons in a B-coherence length size 
region, the impurities can co-self-localize into a self-localized mult i- impurity polaron 
(at least in the parameter region in which the BEC density response to the impurities 
remains small). Further collisions with more impurity atoms grow the self- localized 
polaron impurity numbers. Finally, the relative B BEC density response becomes large, 
the B-BEC density is expelled from the impurity region and the impurity extent grows 
to exceed the B-coherence length. Below, we back up this picture by showing that the 
condition for first order impurity nucleation coincides (or nearly coincides) with the 
condition for zero temperature BEC phase separation. 

5.4. The limit of large co-self-localizing impurity numbers 

As the impurity number A^ increases at fixed (3 coupling, the BEC density response to 
the co-self- localizing impurity polaron grows. For an impurity number A^d, the BEC 
density becomes comparable to the average BEC-density and the self-localized impurity 
structure crosses over to a 'bubble' instead of a sub-coherence length polaron. In that 
regime, the Frohlich-coupling of equation (l22ll . which was derived in the Bogoliubov 
approximation and assumes low depletion, does not accurately describe the impurity- 
boson coupling and the Landau-Pekar treatment breaks down. 

Hence, the large impurity limit, A^ — )■ 00, of /3c further strengthens the connection 
between mult i- impurity polaron formation and BEC phase separation. Keeping only the 
A^^-terms of the E of equation (l52l) , the kinetic energy term drops out and the polaron 
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size is determined from the competition of the BEC B mediated interactions and the 
short-range impurity-impurity interactions, 



E 



(N) 



V 

2x3 



(62) 



M exp ^— ^ j exp ^— u| j du. Introducing y = x/(3 as variable, 



where f{x/(3) 

y = a/^, the energy functional takes the form 
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En 
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(63) 



where A = iV^/ [2/3^?7] > and k = -y/S/vr/?^?]"^. This function has a minimum 
in y and the local minimum has a negative value provided k exceeds a minimal 
value Kc- Numerically, we found Kc = 1.0004. When k > Hc, the value of /i/ 
at the minimum E is negative as fii = 2E/N in the absence of the kinetic energy 
term. Therefore, a large number of bosonic impurities co-self-localize ii k > Hc or 
/3 > (tt/S)^^^ y/rjy/K^ = /3c{N — )■ oo). We assume that the nucleation occurs when a 
number, any number, of impurities can co-self- localize. As Pc{N) decreases with A^, 
{N — )■ oo) represents the minimal coupling value at which the transition takes place. 
Therefore, a transition sets in if 
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/3 > (tt/S)^/^ v^v^ = (16vr)^/^ 

V (iBB 



where we used r] = 8^/2(3'j{aii/aBB)- As a consequence, the nucleation condition takes 
the form 

' an 
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P > 4v/^/€,7 ( ^ ) (65) 

V PoO'%B^ obtain 

O-IldBB 4 \ TTlB J \ mj J 

which is very close to the phase separation condition of equation (!6T]) as = 1.0004. 

What is Nd at which the self-localized polaron structure crosses over to a 
phase separated bubble? The peak BEC density variation for a Gaussian impurity 
wavefunction takes the form 
5pij(x 
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Po 7rV4 \ mj x' 

In the pinprick regime, ^ 1 with weak impurity- impurity interactions, A^(A^ — 1) -C 
Sa/tiTs/ [27/] j so that x = 3y^7i/2/N, we obtain 

1/4 
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Estimating from |5pB (x = 0) /po| = 1, we find 

(69) 



'37r3/4 fmiV^^ 1 



The A^^-scaling of Spb/po may convey the message that muhi- impurity polarons with 
a finite number of impurities is unhkely. However, at larger impurity number, the self- 
localization also sets in at a smaller value of A^. The relative BEC density variation 
near self- localization with /3 ^ /J^ ~ (27r) /N is 



SpB (x = 0) 



^^tt^/^v/^V^t'^'A^'/' (70) 



PO N,l3=fSc(N) 3 

which scales as \/iV and can be satisfied experimentally for TV 7^ 1. 
5.5. Self-localized impurity polaron fluid 

It is interesting to note that if an is replaced by ^ in the phase separation coupling 
constant the expression goes over into the expression for the impurity-boson coupling 
constant /9, 

sf-P. (71) 

At phase separation, g = 1 and /3 = /3c(A^ 00) so that, up to a factor ^J~k^\ 

^ = fi^N ^ 00) . (72) 

From figure m it is clear that when 77 is sufficiently large to make fiJ^N — )■ 00) comparable 
to or greater than 5, (implying an ^ ,^), then /3c(l) < /3c (A^ > !)• What is the 
ground state of an impurity-BEC system with multiple impurities and an impurity-boson 
coupling constant /3 with /3c(l) < /3 < /3c(A^ — >■ C)o)? The conditions suggest that one 
impurity atom self-localizes but that the co-self-localization of two or more impurities 
is prevented by the short-range impurity-impurity repulsion. Each impurity could self- 
localize, but how do the self-localized single impurity polarons arrange themselves? As 
?7 is large, two self-localized polarons repel each other at small distances but attract 
each other at larger distances ranging up to ^. In classical particle systems, these are 
the conditions at which a crystal can form or a liquid of self-determined density. Would 
the self-localized polaron liquid form a self-localized polaron crystal, liquid or BEC? 

This regime is, however, not easily realized experimentally. In the absence of any 
resonances, the cold atom BEC coherence length generally exceeds the scattering length 
by one or two orders of magnitude. The realization of a self-localized single impurity 
polaron liquid that does not nucleate into a phase separated BEC, would require the 
enhancement of both the ajs and the ajj. If one Feshbach resonance enhances aiB to 
increase /3, the magnetic field is 'used up' and an cannot be varied magnetically. We 
suggest that in two-dimensional systems confinement-induced resonances can be used 
to enhance the other scattering length. 
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The confinement of atoms to one or two dimensions by squeezing the trap potential 
in the other dimensions induces resonances in the atom-atom interactions. These occur 
when the extent of the wavefunction in the transverse (squeezed) direction becomes 
comparable to the scattering length in three dimensions [32], [33]. As we expect 
one and two-dimensional BEC or quasi-BEC systems to undergo phase separation as 
well as exhibit self-localization of impurity atoms, we suggest that a combination of a 
magnetically controlled Feshbach resonance and a confinement-induced resonance may 
access self-localized single impurity polaron fiuids. 

6. Conclusions 

The experimental observation of multi-impurity self-localized polarons would provide an 
interesting class of structures that have a subcoherence length size. BEC experiments 
can address fundamental, strong coupling questions: Does polaron self-localization give 
transition like behavior? What is the effective mass of self-localized BEC impurity 
polarons? Can the impurity extent, as well as the polaron size be measured directly 
by species specific dipole potentials? In addition to fundamental many-body research 
the realization of these structures may bring new methods to bear for manipulating 
and probing BEC systems. Self-localized multi-impurity polarons can be transported 
by a laser beam of moving but broad focus. Can the movable self-localized polaron 
carry qubit information? The creation of shock wave analogues in BEC's required 
density perturbations that were smaller than the BEC-coherence length [M], can other 
phenomena be triggered such as superradiance in phonon radiation? The realization 
of a Josephson junction in BEC's described by the usual linear Josephson coupling 
Hamiltonian requires a constriction that is smaller than the BEC-coherence length. 
Can a trapped self-localized impurity provide such constriction? 

If many impurity atoms are present in the trap, the self-localized impurity polarons 
may represent transient structures. The self-localization dynamics, we suggest, describes 
the onset of the phase separation nucleation process. On the other hand, existing cold 
atom techniques can isolate the self-localized impurity droplet with a controlled number 
of impurities. For instance, a species specific optical lattice can contain a well determined 
number of impurity atoms per site. An overlapping BEC does not feel the optical 
lattice but can self-localize the A^-impurities. As self-localization depends on the number 
of co-located impurity atoms, this technique could provide another sensitive method for 
counting atoms in each well. 

In summary, we have shown that N bosonic neutral impurity atoms embedded 
in a dilute gas BEC can self-localize into a collective strongly coupled polaron state. 
The critical impurity-boson interaction strength for self-localization can be significantly 
lower than the corresponding value for a single impurity polaron. We identify the 
self-localized impurity polarons as the initial droplets of the BEC phase separation 
nucleation process. In the limit of large impurity- impurity interactions, there may be a 
strongly interacting regime in which self-localized single impurity polarons form a stable 
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self-localized polaron liquid. We have discussed the experimental challenge of realizing 
these systems in the polaron regime and not as phase separated bubbles. Finally, we 
have speculated on possible applications of this new class of sub-coherence length BEC 
structure. 
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